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The factorization theorem for exclusive processes in perturbative QCD predicts the behavior 
of the pion electromagnetic form factor F{t) at asymptotic spacelike momenta t(= —Q^) < 0. 
We address the question of the onset energy using a suitable mathematical framework of analytic 
continuation, which uses as input the phase of the form factor below the first inelastic threshold, 
known with great precision through the Fermi- Watson theorem from tttt elastic scattering, and 
the modulus measured from threshold up to 3 GeV by the BaBar Collaboration. The method 
leads to almost model-independent upper and lower bounds on the spacelike form factor. Further 
inclusion of the value of the charge radius and the experimental value at —2.45 GeV^ measured 
at JLab considerably increases the strength of the bounds in the region Q^ < lOGeV^, excluding 
the onset of the asymptotic perturbative QCD regime for Q^ < 7GeV^. We also compare the 
bounds with available experimental data and with several theoretical models proposed for the low 
and intermediate spacelike region. 
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I. INTRODUCTION 



The high-energy behavior along the spacelike axis of 
the pion electromagnetic form factor is predicted by fac- 
torization in perturbative QCD [ll-U, which to leading 
order (LO) gives"^ 



-^pertl ^ 



STTf^asili') 



(1) 



where /^r = 130.4 MeV is the pion decay constant 
and as{^i'^) the strong coupling at the renormalization 
scale ^^. Next-to- leading-order (NLO) perturbative cor- 
rections have been calculated in |5l-[ll|. using various 
renormalization schemes and pion distribution ampli- 
tudes (DAs). In particular, the result in the MS- 
renormalization scheme with asymptotic DAs reads [10| 
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l3o A /^' 



14 
-1-M2 



(2) 
where /3o = 11 — "iuf/S is the first coefficient in the 
perturbative expansion of the /3-function, n/ being the 
number of active flavors. The ambiguities that affect the 
perturbative QCD predictions have been investigated in 
many papers, where, in particular, the dependence on 



^ In our convention the form factor F{t) is a function of the squared 
momentum transfer t = q^, and is a real analytic function in 
the complex i-plane cut along t > 4M^. In this notation the 
spacelike axis is defined by t(= —Q^) < 0. 



the renormalization scale fj? and various prescriptions 
for scale setting have been discussed [ifl, [l3 • 

It has been known for a long time that in the case of 
the pion form factor the asymptotic regime sets in quite 
slowly, due to the complexity of soft, nonperturbative 
processes in QCD in the intermediate Q^ region. Many 
nonperturbative approaches have been proposed for the 
study of the form factor, inclu ding QCD sum rules [l3| . 
quark- hadron local duality [l4J17| . extended vector me- 
son dominance fl8'|, light-cone sum rules [19l - l2]| . sum 
rules with nonlocal condensates |22h241 |. and anti-de Sit- 
ter (AdS)/QCD models [H,!!!]. The scale of the onset 
in the presence of Sudakov corrections [27| and large Nc 
Regge approaches [2^] suggests that it can be very large. 
However, constructing a fully valid model to describe the 
form factor at intermediate energies in fundamental QCD 
still remains a major theoretical challenge. 

On the experimental side, several measurements of the 
spacelike form factor at various energies are available (see 
Refs. ,29h39i])j the most precise being the recent results 
of the JLab Collaboration |3^, |30| with the highest mea- 
surement being at Q^ = 2.45 GeV^. The experimental 
determination of F at larger values of Q^ is difficult due 
to the lack of a free pion target and requires the use of 
pion electroproduction from a nucleon target. This stems 
from a virtual photon coupling to a pion in the cloud sur- 
rounding the nucleon. In this regard, there are uncertain- 
ties associated with the off shellness of the struck pion 
and the consequent extrapolation to the physical pion 
mass pole, which leads to uncertainties in the extraction 
of the value of the form factor at t = 0, and possible 
contributions from nonresonant backgrounds. The lack 
of reliable experimental data in the higher Q^ region is a 



major obstacle to confirm or discard the number of theo- 
retical models already available. In particular, it remains 
an open question as to at what value of Q^ do the non- 
perturbative contributions become negligible so that the 
perturbative QCD description of the form factor becomes 
reliable. 

In the recent past, the knowledge of the form factor has 
improved considerably on the timelike axis. The phase 
of the form factor in the elastic region of the unitarity 
cut (below the first inelastic threshold associated with 
the TTtLi state) is equal, by the final state interaction the- 
orem also known as the Fermi- Watson theorem, to the 
P-wave phase shift of the elastic tttt scattering ampli- 
tude, which has been calculated recently with precision 
using Roy equations pOl - l43l | . The modulus has been mea- 
sured from the cross section of e^ e~ — >■ 7r"''7r~ by several 
groups in the past, and more recently to high accuracy 
by BaBar [ii] and KLOE [45", "iS] Collaborations. In par- 
ticular, the high statistics measurement by BaBar yields 
the modulus at high precision up to an energy of 3 GeV. 
Therefore, an interesting possibility is to find constraints 
on the spacelike form factor by the analytic extrapolation 
of the accurate timelike information. 

In the present paper we address precisely this problem. 
Namely, our aim is to perform an analytic continuation 
of the form factor from the timelike to the spacelike re- 
gion using in a most conservative way the available infor- 
mation on the phase and modulus on the unitarity cut, 
and also spacelike information. The main result of the 
method will consist of rather stringent upper and lower 
bounds at values of spacelike momenta, which provide 
a criterion for finding a lower limit for the onset of the 
QCD perturbative behavior. 

We start by briefly discussing, in Sec|lTl several meth- 
ods of analysis based on the analytic properties of the 
form factor investigated in the literature [47- 5q. We 
then present our method, which uses the available input 
in a very conservative way. As we will explain, the price 
to be paid is that we are able to derive only upper and 
lower bounds on the spacelike form factor. The extremal 
problem leading to these bounds is formulated accord- 
ingly 

In Sec lIIII we present the solution of the extremal prob- 
lem formulated in the previous section. We use a math- 
ematical technique applied for the first time in [50| and 
reviewed more recently in [60] , which has been exploited 
also for the investigation of the low-energy shape param- 
eters and the location of the zeros of several hadronic 
form factors f60|-t64|. We emphasize that the solution of 
the problem is exact, so within the adopted assumptions 
the bounds are optimal. 

In Sec II VI we gather all the theoretical and phenomeno- 
logical inputs that go into our calculation, and in SecJVl 
we present the bounds on the form factor in the space- 
like region derived using our formalism. In this section 
we also compare our findings with both experimental and 
theoretical results available, being able, in particular, to 
investigate the onset of the asymptotic regime of pertur- 



bative QCD and to check the validity of various nonper- 
turbative models proposed in the literature. Finally in 
Sec I VII we present a brief discussion on the implication 
of our results and draw our conclusions. 



II. ANALYTIC CONTINUATION 

The standard dispersion relation in terms of the imagi- 
nary part on the form factor along the unitarity cut is the 
simplest way to perform the analytic continuation from 
the timelike to the spacelike region. This method was 



applied for a discussion of perturbative QCD in [47[ , and 
more recently in J48{ . As the imaginary part of the form 
factor is not directly measurable, the method requires 
the knowledge of both the phase and the modulus, and 
is influenced by errors coming from the relatively poor 
determination of the modulus at low energies, and the 
lack of knowledge of the phase above the elastic region. 

Alternative dispersive analyses applied so far are based 
on representations either in terms of the phase, the so- 
called Omncs representation |49l - [5l| . or in terms of the 
modulus [53 . Both approaches require the knowledge of 
zeros of the form factor, and are plagued, respectively, by 
uncertainties related to the unknown phase in the inelas- 
tic region, and the uncertainties of the modulus at low 
energies. Various analytic representations or expansions 
in terms of suitable sets of functions were also used in 
[5J]-[53| for data analysis and the analytic extrapolation 
of the form factor. 

Our approach aims to use as input in a most conserva- 
tive way the precise information available on the unitarity 
cut. First, we consider the relation 



Arg[F(t + ie)\ = 6\{t), AMI < t < U, 



(3) 



where 6l (t) is the phase shift of the P wave of tttt elas- 
tic scattering. We denoted by ti„ the upper limit of the 
elastic region, which can be taken as tjn = (M^ -I- M(^)^ 
since the first important inelastic threshold in the uni- 
tarity relation for the pion form factor is due to the ujtt 
pair. 

Under weak assumptions, the asymptotic behavior ([T|) 
along the spacelike axis implies a decrease like \F{t)\ ~ 
1/t also on the timelike axis.^ Therefore, using the recent 
experimental data on the modulus up to \/t = 3 GeV 
[44| . supplemented with conservative assumptions above 
this energy, we can obtain a rather accurate estimate of 
an integral of modulus squared from tin to infinity. More 
precisely, we assume the following condition. 



1 



dtp{t)\F{t)\^^I, 



(4) 



tin 



^ We quote a rigorous result given in l65ll . which states that, if 
an analytic function F{t) has a bounded phase for t — >■ — 00 and 
t — > 00, the ratio \F{—t)/F{t)\ tends asymptotically to 1, at least 
in an averaged sense. 



where p{t) is a suitable positive-definite weight, for which 
the integral converges, and the number / can be esti- 
mated with sufficient precision. The optimal procedure 
is to vary p{t) over a suitable admissible class and take 
the best result. In principle, a large class of positive 
weights, leading to a convergent integral for \F{t)\ com- 
patible with the asymptotic behavior (jlj of the pion form 
factor, can be adopted. A suitable choice is 



Pb,c(i,Qo) 



t" 



(^ 



QIY 



(5) 



with the parameter Qq > and b, c in the range < 
c — 6 < 2. In most of our calculations we shall choose the 
simpler form 



Pa[t) = -, 



(6) 



with < a < 2. The constraint ^ and the choice of the 
optimal weight will be discussed in more detail later. 

Additional information inside the analyticity domain 
can be implemented exactly. In practice we shall use the 
input 



F{Q) = 1, 



F'{Q) = \{rl 



(7) 



with the charge radius (rj) varied within reasonable lim- 
its, and the values of the form factor at some spacelike 
values 



F(t„) = F„±e, 



U < 0, 



(8) 



where Fn and e„ represent the central value and the ex- 
perimental uncertainty, known from the most precise ex- 
periments [38ll39l|. 

The relations dS]), (g]), d?]) and © define a specified 
class of real analytic functions in the i-plane cut for t > 
AM^. The problem is to derive rigorous upper and lower 
bounds on F{t) in the region i < 0, for functions F{t) 
belonging to this class. This extremal problem will be 
solved exactly in the next section. 



III. SOLUTION OF THE MATHEMATICAL 
PROBLEM 

For solving the problem formulated in the previous sec- 
tion, we apply a standard mathematical method [59l |60| . 
We first define the Omnes function 



0(t) = exp 



^ Jam^ t'{t' - t) 



(9) 



where 6{t) — Sl{t) for t < tin, and is an arbitrary func- 
tion, sufficiently smooth {i.e., Lipschitz continuous) for 
t > tin. As shown in [60|, the results do not depend on 
the choice of the function S{t) for t > tjn. The crucial 
remark is that the function h(t) defined by 



is analytic in the i-plane cut only for t > Un- Further- 
more, equality (g]) implies that h{t) satisfies the condition 



dtp{t)\0{t)\'\h{t)\'=l. 



(11) 



tin 



This relation can be written in a canonical form, if we 
perform the conformal transformation 



m = 



''in V ''in '' 



t 



iin + V ^ii 



(12) 



which maps the complex t-plane cut for t > ijn onto the 
unit disk \z\ < 1 in the z plane defined hy z = z{t), and 
define a function g(z) by 



g{z)^wiz)L,{z)F{i{z))[Oii{z))Y 



(13) 



where t{z) is the inverse of z = z(t), for z{t) as defined in 
(J12p . and w{z) and a;(z) are calculable outer functions, 
i.e., functions analytic and without zeros for \z\ < 1, de- 
fined in terms of their modulus on the boundary, related 
to \/p{t) and |0(t)|, respectively. 

It follows from (fTO)) that the product 
F{i{z))[0{i{z))]^^ appearing in ^ is equal to 
the function h(t{z)) , which is analytic in j^j < 1. 
Therefore, the function g{z) defined in (|13p is analytic 
in \z\ < 1. 

The outer functions corresponding to the weight func- 
tions ^ and ^ can be written in an analytic closed 
form in the z variable as 



Wb^ciz,Qo) = C^VU^) 



l+b-c 



(l-z)V^ {l + ~z{-Ql)r 

(l + z)3/2-c+b(l_22(_g2))c 

(14) 



and, respectively, 

Waiz) = i2y/t~) 



1-a 



(1 



U/2 



(l + z)3/2-a- 



(15) 



For the outer function a; we shall use an integral rep- 
resentation in terms of its modulus on the cut t > tin, 
which can be written as [59, roOl 



uj{z) — exp 



\/tin -i{z) 



\n\0{t')\dt' 



t,„ ^F^t~{t' -t{z)) 



(16) 
Using the definition of the functions w{z) and uj{z), it 
is easy to see that (ITT]) can be written in terms of the 
function g{z) defined in p^ as 



1 r'^'^ 



C^ 



r,iS 



(17) 



We further note that p2p implies that the origin t ~ 
of the t plane is mapped onto the origin z = of the z 
plane. Therefore, from (TlBl it follows that each coefficient 
gk (z R oi the expansion 



F{t) = 0{t)h{t) 



(10) 



9{z) = 50 + giz + g2z^ + 93z^ + 



(18) 



is expressed in terms of the coefBcients of order lower or 
equal to fc, of the Taylor series expansion of the form 
factor at t — 0. Moreover, the values F{tn) of the form 
factor at a set of real points t„ < 0, n = 1, 2, ..., N lead 
to the values 



g{Zn) = w{Zn)uj{Zn)F{tn) [0{tn)\ 



z{tn)- 



(19) 

Then the L^ norm condition (J17p implies the determi- 
nantal inequality (for a proof and older references see 






^1 



{Z1Z2 
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N 
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K 
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{ziZm)'' 
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,2K 



'JV 



>0, (20) 



where 



K-1 



K-1 



(21) 



fe=0 



fc=0 



K > 1 denoting the number of Taylor coefficients from 
(TTSl) included as input. 

The inequality ((20)) leads to rigorous bounds on the 
value of the form factor at one spacelike point using input 
values at other spacelike points. We can implement also 
a number K of derivatives at i = 0, in particular, the nor- 
malization and the charge radius from ^. The deriva- 
tion of the upper and lower bounds amounts to solving 
simple quadratic equations for the quantities g{zn), re- 
lated to the values F{tn) of the form factor through the 
relation P^ . As shown in [GO"], the inequality (|20p holds 
also if the equality sign in ^ is replaced by the less than 
or equal sign. Furthermore, as argued in |60|, the bounds 
depend in a monotonic way on the value of the quantity 
/ appearing in (1171) . becoming weaker when this value is 
increased. 



IV. INPUT 

The phase shift 5\{t) was determined recently with 
high precision from Roy equations applied to the tttt elas- 
tic amplitude in |40l - t43l l. We use as phenomenological 
input the phase parametrized in [43] by 



^ot5\{t)^^^{Ml-A 



2Ml 



APpVt 



Ba + B, 



where k = y^t/A- M^, 



= 1.05 GeV, 



(22) 
Bo = 1.043 ± 



0.011, Bi = 0.19 ± 0.05 and Mp = 773.6 ± 0.9 MeV. We 
assume isospin symmetry and take for Af^r the mass of 
the charged pion. 



The function 5\ obtained from (P^ with the central 
values of the parameters is practically identical with the 
phase shift obtained in [40| from Roy equations for \/t < 
0.8 GeV, and its uncertainty in the whole elastic region 
t < ijn is very small. 

Above tin — (0.917 GeV)^ we use in Q a continuous 
function S(t) , which approaches asymptotically n. As 
shown in [60| , if this function is sufficiently smooth (more 
exactly, Lipschitz continuous), the dependence on S{t) of 
the functions O and uj, defined in (O and ([TO)), respec- 
tively, exactly compensate each other, leading to results 
fully independent of the unknown phase in the inelastic 
region. 

For the calculation of the integral defined in Q we 
have used the BaBar data J44'] from \/t\^ — 0.917 GeV 
up to ^/t = 3 GeV, and have taken a constant value for 
the modulus in the range 3 GeV < \/t < 20 GeV, con- 
tinued with a 1/t decrease above 20 GeV. This model 
is expected to overestimate the true value of the inte- 
gral: indeed, we take up to 20 GeV the modulus equal 
to 0.066, i.e., the central BaBar value at 3 GeV, while 
the perturbative QCD expression ([T]), continued to the 
timelike axis, predicts a much lower modulus, equal to 
0.011 at 3 GeV and to 0.000 16 at 20 GeV. According to 
the above discussion, a larger value of / leads to weaker 
bounds. Therefore, if we use an overestimate of the in- 
tegral, we weaken the bounds which nevertheless remain 
valid. This makes our procedure very robust. 

As we mentioned, in our analysis we shall work with 
weights of the form ([6]). The values of / corresponding 
to several choices of the parameter a are given in Table 
m where the uncertainties are due to the BaBar exper- 
imental errors. As expected, the most sensitive to the 
uncertainty of the high-energy behavior of |i^(t)| is the 
integral / for a = 0, because in this case the weight is 
constant. We checked, for instance, that if we use in the 
integral above 20 GeV the modulus \F{t)\ equal to 0.000 
16 instead of 0.066, we obtain for / a value smaller by 
about 24% for a — 0, while the change in / is of only 3% 
for a = 1/2 , and for a = 1 and a = 2 the values of / 
remain practically unchanged. 

We complete the specification of our input by giving 
the range [5^, [53| 



(rl) 



0.43 ± 0.01 fm^ 



(23) 



adopted for the charge radius, and the spacelike datum 
38, 391 



F(-2.45GeV2) 



0.167 iO.OlOtaoo?- 



(24) 



We used the spacelike datum which gives the most 
stringent constraints. In the present case this is the point 
that is closest to the higher energies of interest. In con- 
trast, in the case of the shape parameters, the space- 
like data closest to the origin provide the most stringent 
constraints [63] ■ Of course, more spacelike data may be 
included as input but, as discussed in [63, when the un- 
certainties are taken into account the improvement is not 
significant. 



TABLE I: Values of the integral / defined in Q corresponding 
to different weight functions defined in ((6]). 
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FIG. 1: Upper and lower bounds obtained with several weight 
functions pa{t) ~ l/t". No spacelike information is used as 
input. The black lines correspond to the central values of the 
input and the grey bands indicate the errors. 



V. RESULTS 

For illustration we show first, in Fig. [TJ the upper and 
lower bounds on the product Q^ F{—Q'^) in the range 
< Q^ < 10 GeV^, obtained with different weight func- 
tions, using as input the relations (jSj, (01) and (O, with 
no information from the spacelike axis. The solid lines 
are the bounds for the central values of the input, and the 
grey bands indicate the uncertainty on the corresponding 
bounds, obtained by adding in quadrature the uncertain- 
ties due to the variation of the phase given in ( 1^ , the 



charge radius (r^) given in (|23p . and the integral / from 
Table H 

The results obtained with different weights are rather 
similar, except for the weight P2(i), which gives much 
weaker bounds at high Q^ . This feature is actually ex- 
pected: since this weight decreases rapidly, the class of 
functions satisfying Q is less constrained at high ener- 
gies, even functions increasing asymptotically being ac- 
cepted. We note that in the so-called unitarity bounds 
approach, applied to the pion form factor in [5^, [6j] , an 
inequality similar to Q is derived from a dispersion rela- 
tion for a physical observable, like a polarization function 
calculated in perturbative QCD or the hadronic contri- 
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FIG. 2: Effect of the additional spacelike input (|24p on the 
bounds obtained with the weight pi/2(t). Central values of 
all the input parameters are used. 
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FIG. 3: Upper and lower bounds obtained with various 
weights Paify = l/i", for the central values of the input pa- 
rameters. The inset shows the bounds up to 1 GeV^. 



bution to the muon anomaly, using in addition unitarity 
for the corresponding spectral function. In this approach 
the effective weight p(i) is fixed and cannot be chosen at 
will, decreasing actually like Xjt^ in all the cases inves- 
tigated. Therefore, as discussed in [63j] , the approach of 
unitarity bounds, which is useful when no information on 
the modulus is available, would give weak bounds on the 
pion form factor at higher energies on the spacelike axis. 
The recent high statistic measurements of the modulus 
up to rather high energies on the timelike axis [4^ allow 
us to choose an optimal weight, leading to a remarkable 
improvement of the bounds on the spacelike axis. 

In Fig. [5]we demonstrate the effect of an additional in- 
put from the spacelike datum (j24p . using for illustration 
the weight pi/2(0- The solid lines show the bounds ob- 
tained as before without any spacelike information, while 



the dashed hnes are obtained by imposing the central 
value from ([M)) . The inclusion of this additional infor- 
mation narrows considerably the allowed domain situated 
between the upper and the lower bounds. 

The comparison of various weights is seen in Fig. |31 
where the bounds are again obtained with the central in- 
put values, including the spacelike datum from ([M)) . At 
low Q^ the bounds obtained with different weights are 
almost identical, and are very tight. At higher Q^ the 
small constraining effect of the weight p2 (0 i^ again visi- 
ble, while the other weights give comparable results. By 
inspecting the curves, we conclude that the best results 
are obtained with the weight pi/2{t), and we shall adopt 
this weight in what follows. In fact, it gives bounds al- 
most identical with the weight pi(i) = 1. In addition, it 
has the advantage of being less sensitive to the unknown 
high-energy behavior of the form factor, as discussed in 
the previous section in connection with the calculation of 
the quantities / given in Table ID 

In Fig. |4]we illustrate, for the special weight pi/2{t), 
the effect of the uncertainties of the input on the resulting 
bounds. Since we are interested in the most conservative 
results, i.e., the largest allowed domain defined by the 
upper and the lower bounds, it is enough to exhibit the 
upward shift of the upper bound, and the downward shift 
of the lower bound, produced by the uncertainties of the 
input. By including these uncertainties, the allowed do- 
main obtained with the central input, shown as the inner 
white region in FigUl is enlarged covering the large grey 
domains. We recall that in each point the error is ob- 
tained by adding quadratically the errors produced by 
the variation of the phase (|22|) . the charge radius (l23l) . 
the integral / for a = 1/2 from Table HI and the experi- 
mental value given in ([M]) . 

It turns out that the greatest contribution to the size of 
the grey domain is the experimental uncertainty of the 
spacelike value (PH, which decreases the lower bound 
by about 30% for values of Q^ around 7GeV^, for in- 
stance. This shows that more accurate experimental data 
or lattice calculations at some points on the spacelike 
axis will lead to more stringent constraints in this re- 
gion, which, as shown below, is of interest for the onset 
of the perturbative QCD regime. At larger values of Q^, 
the effect of the uncertainty in the charge radius starts to 
increase, becoming of comparable value, of about 30%, 
near Q^ — lOGeV^. The uncertainty of the integral / 
has an effect of less than 8% in the considered range (as 
expected, it has a noticeable effect only at higher mo- 
menta). 

In Figs. [S] and [51 we compare the constraints corre- 
sponding to the weight p^ /-jit) with some of the data 
available from experiments [22-131, [sHH, [H, ill. We 
find that at low Q^ most of the low-energy data of Amen- 
dolia [35| are consistent with the narrow allowed band 
predicted by our analysis. The most recent data from 
|38l |39| are well accommodated within our band, which 
is not surprising as we used one of these points as input. 
There are however some inconsistencies between the al- 



lowed domain derived here and the data of Bebek et al 
[30, l31| , Ackermann et al [32] , and Brauel et al [3J| , in 
spite of their rather large errors. 




FIG. 4: Upper and lower bounds derived with the weight 
Pi/2(i) by including the charge radius (I23|) and the spacelike 
datum (|24p . Inner white region: allowed domain for the cen- 
tral values of the input. Grey bands: enlarged allowed domain 
with inclusion of the errors of the input. 
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FIG. 5: Allowed domain obtained with the weight pi/2(i) 
compared with several sets of experimental data. 



Finally, in Fig. [7] we compare the allowed domain ob- 
tained in this work with the predictions of perturbative 
QCD and several nonperturbative models. We show first 
the LO expression ^, in which we have taken the scale 
fj,^ — Q^ and used the running coupling to one loop 



a.(Q') 



47r 



/3oln(gVA2)- 



(25) 



We have taken Uf = 3 active fiavors and A = 0.214 GeV, 
which in Eq. (p5|) gives as{M^) = 0.33, the average of 
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FIG. 7: Comparison of the bounds with perturbative QCD 
and several nonperturbative models. 



the vario us p redictions from hadronic r decays obtained 
recently [66l - l72l |. As seen from Fig. [71 the central LO 
curve, obtained with this choice of the coupling, stays 
below the lower bound derived in the present work up to 
Q"^ — 7 GeV^. The variations set by the error of about 
±0.02 in as(M2) quoted in [H-lzl do not modify this 
conclusion. 

As discussed in [10, [l2l the perturbative prediction 
to NLO is sensitive to the choice of the renormalization 
scale, and also of the factorization scale in the case when 
pion DAs different from the asymptotic ones are used 
in the calculation. Several prescriptions for scale setting 
were adopted, but there is no general consensus on the 
issue. In particular, in the Brodsky-Lepage-Mackenzie 
procedure the renormalization scale is chosen such that 
the coefficient of /3o in the NLO correction ([2]) vanishes 

For illustration, in Fig. [7] we show the sum of the 
LO and NLO terms ^ and ([2]), obtained with the scale 



Q^ and the one loop coupling (1^ . As in previous 
references [3, [lO, 120, S, |43] , we have taken the number 
of active flavors nj = 3. This curve is compatible with 
our bounds enlarged by errors only for Q^ > 6 GeV . 

When Q^ is varied over a large range, one must of 
course take into account the transition from Uf to Uf + 1 
active flavors in the /3-function coefficients when the 
threshold of a new flavor is crossed. However, it turns 
out that in the region of interest to us, Q^ < 10 GeV^, 
the effect of the higher mass quarks is small. First, we 
recall that the definition of the flavor matching scale is 
subject to a certain uncertainty. If the thresholds are 
set at the pole masses of the charm and bottom quarks, 
nf should be raised to 4 and 5, respectively, when the 
thresholds of 1.65 GeV and 4.75 GeV are crossed. How- 
ever, the change in /3o is small, from /3o = 9 to /3o = 8.33 
for the charm threshold relevant in the region consid- 
ered. Moreover, the continuity of the coupling (j25p must 
be imposed by a change of the parameter A above the 
threshold. Therefore, the modification of the form factor 
below Q^ = 10 GeV is expected to be quite small and 
the comparison with the bounds derived in the present 
paper remains unchanged. We note also that for the 
choice of the quark-flavor matching scales at 2mq as in 
6|, or at ^* = 3.729 GeV and ^^ = 10.588 GeV as in 
67| , the thresholds are not crossed and the number of the 
active flavors remains equal to three in the whole region 
shown in Fig. [T] 

In Fig. [7]we show also several nonperturbative models 
proposed in the literature for th e sp acelike form factor at 
intermediate^region [H, |l3, [13, [23^26} 

In Ref. [20], the authors applied light-cone QCD sum 
rules and parametrized with a simple expression the non- 
perturbative correction, to be added to the LO-I-NLO 
perturbative prediction in the region 1 < Q^ < 15 GeV^. 
In Fig. [7] we show the sum of the soft correction and 
the perturbative QCD prediction to NLO, evaluated at 
a scale n'^ = .5 Q'^ + M"^ with ]VP = 1 GeV^ as argued 
in [20] ■ The model is quite compatible with our bounds, 
the corresponding curve being inside the small white in- 
ner domain for Q^ > 6 GeV^. 

The model based on local duality [la] is also consistent 
with the allowed domain derived here for Q^ > 1 GeV^. 
We mention that this model, proposed in [ij, was re- 
cently developed by several authors 17[ . The other mod- 
els shown in Fig. [7]arc consistent with the bounds derived 
by us at low Q^, but are at the upper limit of the allowed 
domain at higher Q^. The agreement is somewhat bet- 
ter for the model discussed in [23], which is a LO-I-NLO 
perturbative calculation using nonasymptotic pion DAs 
evolved to NLO, with a modification of the QCD cou- 
pling by the so-called analytic perturbation theory. The 
AdS/QCD model considered in [26] is in fact a simple 
dipole interpolation, which is valid at low energies but 
seems to overestimate the form factor at larger momenta. 
The same remark holds for the models discussed in |24l ] 
and [25|, based on QCD sum rules with nonlocal con- 
densates, and the chiral limit of the hard- wall AdS/QCD 
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approach, respectively. 



VI. DISCUSSION AND CONCLUSIONS 

In this paper we derived upper and lower bounds on the 
pion electromagnetic form factor along the spacelike axis, 
by exploiting in a conservative way the precise recent 
information on the phase and modulus on complementary 
regions of the timelike axis. 

More exactly, we have applied a method of analytic 
continuation which uses as input the phase only in the 
elastic region t < tin, where it is known with high accu- 
racy from the dispersive theory of tttt scattering via the 
Fermi- Watson theorem. It can be shown rigorously [60] 
that the results are independent of the unknown phase 
of the form factor above the first inelastic threshold ti^- 
As for the modulus, we have used the BaBar data in the 
range between the first inelastic threshold and 3 GeV 
p4| . and very conservative assumptions above 3 GeV. 
The results are not very sensitive to these assumptions, 
since we include the information on the modulus through 
the L^ integral condition ^, instead of imposing this 
knowledge pointwise, at each t. Obviously, the detailed 
behavior of \F{t)\ is averaged in the integral, and the 
high-energy contribution can be suppressed by a suitable 
choice of the weight p{t). In our calculations we have 
choosen a suitable weight, which on the one hand gives 
sufficiently tight bounds, and on the other hand ensures a 
small sensitivity to the high energy part. We have found 
that these constraints are met by a weight of the simple 
form ^, with the choice a = 1/2. 

We mention that, for completeness, we have investi- 
gated also the more general class of weight functions of 
the form ([5|) and found, for instance, that with the choices 
b — 1/2, c— I and Qq in the range 3 — 6 GeV^ the bounds 
are slightly better and have smaller uncertainties than 
those shown in Fig. |4l For the conclusions formulated in 
this paper the improvement is not essential. However, a 
more systematic optimization with respect to the weight 
is of interest for further studies, especially when more 
accurate input from the spacelike axis will be available. 

Besides the use of the modulus above ti„ only in an 
averaged way, as in the condition Q, we note that our 
method does not exploit completely the present knowl- 
edge of the modulus, since we do not use the data on 
\F{t)\ below tin- While it may seem that this additional 
information can improve in a significant way the bounds, 
it turns out that it is not so. Indeed, the knowledge of 
the phase below tjn has, in the case of the pion form 
factor, a considerable constraining power on the modu- 
lus in the same region. The reason is the strong peak 
produced by the p resonance. The Omnes function 0{t) 
defined in ^ reproduces well this impressive behavior of 
the modulus, and this holds for a large class of choices of 



the arbitrary phase above the elastic region. The effect 
of this unknown phase is actually completely removed in 
our formalism by the information on the modulus, which 
constrains the auxiliary function h(t) appearing in the 
expression ^TU\i of F{t). Therefore, imposing additional 
data on \F{t)\ below 0.9 GeV is not expected to improve 
considerably the bounds, if the input values of the phase 
and modulus in this region are consistent. Of course, 
one can reverse the argument and find constraints on the 
modulus on the timelike axis below tin, in order to test 
the consistency. This analysis will be carried out in the 
future. 

In the present work we have investigated the consis- 
tency of the timelike phase and modulus with the experi- 
mental data and theoretical models on the spacelike axis. 
The slow onset of the asymptotic regime predicted by 
factorization and perturbative QCD for the pion electro- 
magnetic form factor has been known for a long time. A 
lot of theoretical work has been done by several groups on 
the determination of the pion form factor in the spacelike 
intermediate energy region, where the soft, nonperturba- 
tive contributions are expected to be important. How- 
ever, definite conclusions on the validity of the theoretical 
models and the precise energy at which the asymptotic 
regime may be considered valid were not possible, due 
to the lack of accurate experimental data in the relevant 
region. 

The bounds derived in the present paper are almost 
model- independent and very robust, allowing us to make 
some definite statements. From Fig. [7] it is possible to 
say with great confidence that perturbative QCD to LO 
is excluded for Q^ < 7 GeV^, and perturbative QCD 
to NLO is excluded for Q^ < 6 GeV^, respectively. If 
we restrict to the inner white allowed domain obtained 
with the central values of the input, the exclusion regions 
become Q^ < 9 GeV^ and Q^ < 8 GeV^, respectively. 
Among the theoretical models, the light-cone QCD sum 
rules [20| and the local quark-hadron duality model [la] 
are consistent with the allowed domain derived here for 
a large energy interval, while the remaining models are 
consistent with the bounds at low energies, but seem to 
predict too high values at higher Q^. 

To increase the strength of the predictions, a reduction 
of the grey bands produced by the uncertainties of the 
input is desirable. As we mentioned, the biggest contri- 
bution is brought by the experimental errors of the input 
spacelike datum ([M| . Therefore, more accurate data at 
a few spacelike points, particularly at larger values of Q^ 
will be very important for increasing the predictive power 
of the formalism developed and applied in the present pa- 
per. 
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